Quantum phase transitions in attractive extended Bose-Hubbard Model with 

three-body constraint 
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The effect of nearest-neighbor repulsion on the ground-state phase diagrams of three-body con- 
strained attractive Bose lattice gases is explored numerically. When the repulsion is turned on, in 
addition to the uniform Mott insulating state and two superfiuid phases (the atomic and the dimer 
superfluids) , a dimer checkerboard solid state appears at unit filling, where boson pairs form a solid 
with checkerboard structure. We find also that the first-order transitions between the uniform Mott 
insulating state and the atomic superfiuid state can be turned into the continuous ones as the re- 
pulsion is increased. Moreover, the stability regions of the dimer superfiuid phase can be extended 
to modest values of the hopping parameter by tuning the strength of the repulsion. Our conclusions 
hence shed light on the search of the dimer superfiuid phase in real ultracold Bose gases in optical 
lattices. 



Spectacular progress both in theories and experiments 
has recently been made on ultracold atomic and molecu- 
lar gases in optical lattices. Owing to the remarkable con- 
trol over physical parameters, ultracold gases offer oppor- 
tunities to simulate the physics of strongly correlated sys- 
tems in regimes which are not easily accessible to solid- 
state materials. As a result, they provide very clean and 
tunable systems in the search for exotic quantum phases 
and in probing quantum critical behaviors around these 
phases^ For instance, successful experimental realization 
of the superfluid-Mott transition for ultracold bosons in 
an optical lattice^ has paved the way for studying other 
strongly correlated phases in various lattice models. 

Very recently, it was suggested that intriguing quan- 
tum critical behaviors could occur in attractive bosonic 
lattice gases with three-body on-site constraint j2» The 
on-site constraint can arise naturally due to large three- 
body loss processes^ and it stabilizes the attractive 
bosonic systems against collapse. Such three-body con- 
strained systems can be realized also in Mott insulat- 
ing states of ultracold spin-one atoms at unit filling^ As 
found by the authors in Ref. 3j, a dimer superfiuid (DSF) 
phase consisting of the condensation of boson pairs can be 
realized under sufficiently strong attraction. According 
to their analysis, the transitions between the DSF phase 
and the conventional atomic superfiuid (ASF) state are 
proposed to be of Ising-like at unit filling and driven first- 
order by fluctuations via the Coleman- Weinberg mech- 
anism® at other fractional fillings. Later investigation 
focuses on the nature of the superfluid-insulator transi- 
tions^ It is shown that, while the Mott-insulator (MI) to 
DSF transitions are always of second order, the continu- 
ous MI- ASF transitions can be preempted by first-order 
ones and interesting tricritical points can thus appear on 
the MI- ASF phase boundaries. The conclusions obtained 
in Refs. and [3 are partly supported by a recent numer- 
ical study employing stochastic series expansion (SSE) 
quantum Monte Carlo (QMC) method implemented with 
a generalized directed loop algorithm.^ In particular, the 
existence of a tricritical point along the saturation tran- 
sition line is verified. However, the nature of MI-DSF 



transitions is not examined in their QMC work. 

In the present work, the effect of the nearest-neighbor 
mutual repulsion on the ground-state phase diagrams 
of three-body constrained attractive lattice bosons is 
investigated by means of exact diagonalizations (ED). 
The systems under consideration are described by the 
extended Bose-Hubbard model with a three-body con- 
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Here, <2j(aj) is the bosonic annihilation (creation) oper- 
ator at site i, t is the nearest-neighbor hopping integral, 
U < is the on-site two-body attraction (\U\ = 1 as the 
energy unit), and u is the chemical potential. V > 
denotes the nearest-neighbor repulsion, which can come 
from the dipole-dipole interactions of the dipolar bosons 
polarized perpendicularly to the lattice plane by truncat- 
ing it at the nearest-neighbor distance. For the discus- 
sions of possible experimental realizations on the present 
model, we refer to Ref. lfj and references therein. The 
averaged particle density is denoted by n and periodic 
boundary conditions are assumed. Our main results are 
summarized in Fig. [TJ When V ^ 0, in addition to the 
MI and two superfiuid phases mentioned above, a dimer 
checkerboard solid (DCS) state consisting of the checker- 
board arrangement of boson pairs emerges and occupies 
the middle part of the phase diagram in the small-i limit. 
The DSF states now appear only in between the uniform 
MI and the DCS states. Moreover, finite repulsion has 
interesting effects on the Ml-superfluid transitions also. 
Because a modest nearest-neighbor repulsion can avoid 
cluster formation and then will suppress phase separa- 
tion, we observe that the segments of the first-order MI- 
ASF transitions on phase boundaries of either the n = 
or the n = 2 MI states (denoted by thick blue dashed 
lines in Fig.[I| can shrink to zero as V is increased. Same 
conclusion has been reached in other context M- Besides, 
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FIG. 1: (Color online) Ground-state phase diagram obtained 
by ED with (a) V = 0, (b)(c) V = 1/2, and (d)(e) V = 1. The 
thin red solid lines indicate the continuous MI-DSF transi- 
tions, and the thick blue solid (dashed) lines show the contin- 
uous (first-order) MI- ASF transitions. The phase boundaries 
of the DCS state are denoted by the green solid circles (lines 
are guides to the eye). The error bar of the determined phase 
boundaries is smaller than the width of the lines or the size 
of the symbols. The schematic phase boundaries between the 
ASF and the DSF phases are added as the thin dotted lines 
for clarity. 



our findings show that the phase boundaries of the MI- 
DSF transitions (denoted by thin red solid lines in Fig. [T]) 
and therefore the stability region of the DSF phase in 
the low-density limit can be extended to modest hopping 
parameters t upon tuning V , such that experimental ex- 
ploration of this interesting state becomes more feasible. 
Our work hence provides a useful guide to the experimen- 
tal search of the DSF phase and the associated quantum 
phase transitions in ultracold Bose gases in optical lat- 
tices. 

The details of our analysis are explained below. Ac- 
cording to the discussions in Ref. [H,, one can determine 
the phase boundaries of either the n — MI state (i.e., 
the empty state) or the n = 2 MI state (i.e., the com- 
pletely filled state) for a given t as follows. For the tran- 



sitions out of the n — MI state, we calculate the ex- 
citation energies (relative to the n = state) £7(1) — fx, 
£7(2) — 2[i, and £7(4) — 4/x of the lowest states within the 
subspaccs of fixed total particle numbers 7V p = 1,2, and 
4, respectively. Here E(N p ) denotes the lowest excita- 
tion energy of £/ebh within the fixed N p subspace. From 
these quantities, the two- and the four-particle binding 
energies, A 2p = £7(2) - 2£7(1) and A 4p = £7(4) - 2£7(2), 
can be obtained. For a given value of the hopping pa- 
rameter t, if both A 2p and A 4p are positive, there ex- 
ists no bound states and the energy gap of the single- 
particle state closes first upon increasing the chemical 
potential fj,. This leads to a continuous MI-ASF transi- 
tion at Mc,omi-asf = E(1)M When A 2p < but A 4p > 0, 
instead, a bound state of two particles appears and its 
gap closes first. Thus a continuous MI-DSF transition 
will occur at /i c ,oMi-DSF = £7(2)/2. Aside from these two 
possibilities, the condition of A 4p < gives a precur- 
sor of instability of boson pairs towards cluster forma- 
tion. That is, phase separation emerges and a first-order 
transition will be observed in varying /i. In the present 
ED analysis, we estimate the first-order transition point 
by £7(4)/4, which provides an upper bound of the exact 
value. Similar discussions apply to the transitions out of 
the n = 2 MI state, where the holes created from the 
n = 2 state take the role played by the particles dis- 
cussed previously. We remind that the first-order transi- 
tion points estimated by the four-hole excitation energies 
will be instead lower bounds of the exact values along the 
saturation transition line. 

Here distinct Ml-superfluid transitions are determined 
by the above method for systems of iV s = 10 x 10 sites. 
The analysis for the V — 1/2 case is described below 
for illustration. The results of various binding energies 
as functions of hopping parameter t for this V are pre- 
sented in Fig. [2] As seen from Fig.[2jb), there is a finite 
region of t within which A 4p < 0. It implies that, on 
the phase boundary of the n = state, the continuous 
MI-DSF (for A 2p < but A 4p > 0) and continuous MI- 
ASF (for both A 2p , A 4p > 0) transitions are separated 
by first-order MI-ASF transitions for 0.15 < t < 0.26. 
On the contrary, the two continuous Ml-superfluid tran- 
sitions on the phase boundary of the n = 2 state should 
meet directly at t « 0.068. The transition points /i c 
for a given t can be determined by using the excitation 
energies E(N P ) as explained above.— According to our 
ED calculations, for V — 1/2, the continuous n = 2 
MI-DSF transitions end at (i E , ^ E ) ~ (0.068, 3.55), while 
(£e, /he) — (0.15, —0.64) for that of the continuous n = 
MI-DSF transition line. There exists also a tricritical 
point (i T ,MT) - (0.26,-1.04) on the n = MI-ASF 
transition line separating the continuous from the first- 
order ones. 

We now turn to the discussions on the phase bound- 
aries of the DCS state. Due to the effect of the nearest- 
neighbor repulsion V, boson pairs can form a solid with 
checkerboard structure at unit filling n — N p /N s = 1 
and lead to the DCS state. In the zero-hopping limit, 
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FIG. 2: (Color online) (a) Two- and four-hole binding ener- 
gies, A2h and A4h, around the n — 2 state and (b) two- and 
four-particle binding energies, A2 P and A4 P , around the n — 
state as functions of t for V = 1/2 and N s = 10 x 10. The 
dashed lines separate the transitions of different characters. 
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FIG. 3: (Color online) Total particle numbers N p as functions 
of chemical potential /i for (a) t = 0.1 and (b) t = 0.3 with 
V = 1/2 and system sizes N a — 4 x 4. The insets illustrates 
the complete changes in particle density n from the n — to 
the n = 2 states as /x increases. 



this DCS state can be stabilized when —1/2 < fi < 
— 1/2 + 2zV with the coordination number z — 4 for 
square lattices. The DCS state can melt into the DSF or 
the ASF states under increasing hopping, and its stabil- 
ity region in [i is expected to reduce to zero as t increases. 
This picture is supported by our numerical calculations 
as seen in Fig. [TJ Due to the limitation in numerics, the 
phase boundaries of the DCS state are estimated from 
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FIG. 4: (Color online) Static structure factor S(k) of the 
n = 1 states for various i's with V = 1/2 and N s — 4 x 4. 



the n = 1 plateaus in the [i-n plots for different i's with 
systems sizes of N s = 4 x 4. For illustration, total par- 
ticle numbers N p and the particle density n as functions 
of chemical potential fi for two different values of t's with 
V = 1/2 are presented in Fig. [31 The large plateaus at 
n = 1 in the fi-n plots clearly indicate the presence of 
the DCS state. Two ends of the plateau give the melting 
transition points in fi for a given t, as depicted in Fig. [T] 
Moreover, distinct characters in various transitions can 
be revealed by the way in which the total particle number 
iVp changes upon varying the chemical potential fir^ In 
the conventional ASF phase, the number of particles will 
increase by 1 when the chemical potential is increased. 
However, in the DSF phase, only pairs of bosons appear 
in the system due to the presence of a pairing gap. Thus 
adding a single boson is forbidden and the jumps in the 
particle number by 2 will be observed as [i is varied. Be- 
sides forming pairs, bosons can become unstable towards 
cluster formation such that the particle number jumps 
by a finite amount in the course of tuning fi. This corre- 
sponds to a first-order transition under the change in the 
chemical potential. As shown in Fig. |3j for t = 0.1, the 
system evolves across a continuous MI-DSF transition 
from the n = state to the DSF state, and then follows 
a first-order transition to the DCS state. For larger t 
(say, t = 0.3), the state after the continuous transition 
from the n = MI state becomes the ASF one instead. 
Adding more bosons by further increasing /i, the system 
can again follows a first-order transition to the DCS state. 

To provide support on the nature of the DCS states 
within the n = 1 plateaus in the fi-n curves, the static 



structure factors 5(k) = (l/iV s 2 )< V. 



l ) of the 



n = 1 states for several hopping parameters i's with 
V = 1/2 and system size N s = 4 x 4 are shown in 
Fig. |U It is found that the static structure factors do 
have peaks at the wave vector k = (jr, tt) and thus signal 
the checkerboard pattern of the boson pairs. As observed 
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FIG. 5: (Color online) Root-mean-square separation r rms of 
the (a) particle (b) hole pairs on N e = 10 X 10 square lat- 
tices as functions of t for various repulsions V . Insets: r rms 
and maximal extent i max in hopping integral for the MI-DSF 
transitions as functions of V . 



from Fig. [4] when t is increased, the peak value S(ir, tt) of 
the structure factor will decrease from its classical value 
S(tt, 7t) = 1 at t = 0. This indicates the quantum melting 
of the DCS state into the DSF or the ASF states under 
increasing hopping. 

Some further microscopic details can be uncovered 
by the ED calculations. Similar to the counterpart of 
fermion pairing, a smooth crossover from the on-site pairs 
to the loosely bound pairs may occur also in the present 
attractive boson systems as the hopping parameter t in- 
creases. In Fig. [SJ evidences supporting this expecta- 
tion are presented. Here the coherence length of the 
boson pairs is estimated by the root-mean-square sep- 



aration r rms = (r 2 ) ) 15 ' 16 which is evaluated under the 
ground state of a single pair within the two-particle or the 
two-hole subspaces. Our results show rapid but smooth 
crossovers from the tightly bound molecules to the loosely 
bound pairs, as long as the DSF phase can be stabilized 
up to modest values of t (say, the V = 1 case for the 
particle pairs and the V = 1/2 case for the hole pairs). 
The general dependence on the repulsion V of the max- 
imal extent i max in hopping integral for the stable DSF 
state in the low-density limit and the corresponding r rms 
is shown in the insets of Fig. [5] We find that larger i max 's 
in general lead to longer r rms 's. The dependence of £ max 
and r rms on V is found to be nonmonotonic, and their 
functional forms shows asymmetry between the cases of 
the particle and the hole pairs. Our conclusions presented 
in Figs. Q] and [5] should be of help in determining optimal 
experimental settings in the search of the DSF phase in 
ultracold Bose gases in optical lattices. 

In summary, the ground-state phase diagrams of the 
three-body constrained extended Bose-Hubbard model 
for various repulsions are investigated. Large plateaus 
at n = 1 in the fi-n curves which show the DCS states 
are observed. Finite repulsions modify the Ml-superfluid 
transitions also. We find that the repulsion V can change 
the first-order MI-ASF transitions into the continuous 
ones and the stability regions of the DSF phase can be 
tuned by V. Therefore, carefully adjusting system pa- 
rameters into the suggested parameter regime are neces- 
sary to find experimentally the interesting DSF phase in 
real ultracold Bose gas in optical lattices. 
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